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THE EFFECT OF ABLATION INJECTION ON RADIATIVE 
AND CONVECTIVE HEATING 

James N. Moss 

NASA Langley Research Center, Hampton, Va. , 23365 

ABSTRACT 

A viscous shock- layer analysis for calculating high energy 
equilibrium flow fields about blunt axi symmetric bodies is 
applied to the problem of massive ablation injection with radia- 
tion transport. A nongray radiation model is used that accounts 
for both line and continuum radiation. The solution method is 
direct and provides both stagnation and downstream solutions. 
Results for shock heated air show that phenolic-nylon injection 
is substantially more effective in reducing the wall radiant flux 
than air injection. Also, for large included body angles, the 
wall radiative flux and the coupled phenolic- nylon injection rate 
do not continue to decrease with increasing distance downstream. 

INTRODUCTION 

The flow field about planetary entry vehicles is very com- 
plex because of the many interactions that occur between ablation 
injection, radiation transfer, and viscous and diffusion effects. 
In addition, no experimental data exist which fully duplicate the 
thermodynamic environment for planetary entry regimes or how the 
environment interacts with the thermal protection material. Con- 
sequently, rigorous numerical analysis of high energy flow fields 
is essential] essential to the understanding of the basic mecha- 
nisms involved in energy transport and to the design of adequate 
thermal protection systems. 

Almost all existing radiative viscous flow analyses are 
limited to stagnation solutions. However, Sutton (Ref. [l]) and 
Chou (Ref. [2]) have recently presented downstream solutions for 
viscous radiating flows. Sutton’s analysis is a direct method 
and is applicable for zero to moderate mass injection while 
Chou's analysis accounts for massive blowing but is an indirect 
method. 

A rigorous viscous shock- layer analysis for calculating high 
energy equilibrium flow fields about blunt axisymmetric bodies at 
zero degree angle of attack is presented herein. The present 
analysis, based on the analysis developed by Davis (Ref. [ 3 ]) for 
a perfect gas, includes mass injection, radiation transfer, dif- 
fusion, and viscous effects. The analysis is direct and provides 
both stagnation and downstream solution for all levels of blow- 
ing. Radiation transfer is calculated with a nongray radiation 
model that accounts for line and continuum radiation. Results 
are presented for earth entry conditions to demonstrate the 
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essential features of this analysis. Particular emphasis is 
given to the effects of mass injection on the stagnation and 
downstream flow. Results for both specified and coupled mass 
injection are presented. The injectants considered are the spe- 
cies for "both air and a phenolic- nylon ablator. This is the 
first direct analysis that provides nonsimilar downstream solu- 
tions for radiative viscous flows for all levels of blowing. 


ANALYSIS 

The viscous shock- layer equations for an equilibrium multi- 
component gas mixture are developed in References [4-] and [53 ♦ 
The same viscous shock- layer equations are solved in the present 
study with the addition of radiation transport. The nondimen- 
sional equations in a body oriented coordinate system (Fig. 1) 



Figure 1. Coordinate system. 
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Elemental continuity: 
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and Q, is the divergence of the radiative flux (Q*Rf|/p*u*3 ) . 
superscript * denotes dimensional quantities while j has a 
value of one for an axisymmetric body and zero for a two- 
dimensional body. The bar denotes mixture values while the sub- 
scripts i and Z denote species and elements, respectively. The 
quantity h is the static enthalpy (h*/p*u*2) # P the density 
£p*/p + ), p the pressure (p*/p*u*), P t he viscosity 
C i tffe elemental mass fraction, Ci the species mass fraction, 

M* the molecular weight, 5i i the number of atoms of element l 
in species i, N the total number of species, R* the universal 
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gas constant, Cp the specific heat, and Np r , N^g, and 
Prandtl, Lewis, and Schmidt numbers, respectively. The reference 
viscosity is the viscosity evaluated at the reference temperature 

U w 2 / C p,°°’ 

The boundary conditions at the shock are calculated by using 
the Rankine-Hugoniot relations. At the wall, the no slip and no 
temperature jump boundary conditions are used; consequently, 
u w = 0. The wall temperature and mass injection rate are either 
specified or calculated. For the calculated conditions, the mass 
injection rate is determined from an energy balance at the flow- 
field- ablator interface. The coupled mass injection rate and 
surface temperature are calculated by iterating the solution of 
the governing flow- field equations and the boundary conditions. 

The total heat transferred to the surface boundary is given, 
in nondimen sional form, as 


q T 
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The quantity K is the mixture frozen thermal conductivity 
(K*/V*ef Cp ot ), Jp the diffusion mass flux of species i 
(j*Ej/h^ef^°° and the mass injection rate [ ( p v*) w /(p*u*) w ] . 

The sub scripts w and - denote quantities at the surface boundary 
and adjacent but below the surface boundary, respectively. The radi- 
ative flux, q r = qj/p*u*3^ and the divergence of the radiative flux, Q, 
are calculated with hhe radiation transport model LRAD - 3, as pre- 
sented in Reference [6] . The radiation transport model is an extended 
version of a coupled line and cont inuum model that was originally 
developed by Wilson (Ref. [7] ) . The LRAD - 3 radiation model accounts 
for the effects of nongray self- absorption, radiative cooling, and the 
contribution of atomic line radiation. Twelve continuum frequency 
bands and nine line bands are used. The species considered for deter- 
mining the radiation transport are: C, H, N, and 0 for both con- 
tinuum and lines and CO, C£> C3, O2, H2, C2H, and e“ for con- 
tinuum. Details of the radiation transport analysis are given in 
References [6] and [7] and, therefore, will not be reproduced here. 

The equilibrium composition is determined by a free energy 
minimization calculation as developed in Reference [8], Thermo- 
dynamic and transport properties are calculated for each species. 

In this study, the binary diffusion coefficient is set equal to 
the diffusion coefficient for atomic carbon diffusing into atomic 
nitrogen. 

The method of solving the governing equations and boundary 
conditions is discussed in detail in Reference [5]; therefore, 
only a brief overview of the solution procedure is presented. 

The governing equations are written in finite-difference form by 
using Taylor's series expansions. A variable grid spacing is 
used in both the tangential and normal directions to the surface 
so that the grid spacing can be made small in the region of large 
gradients. The order of the truncation terms neglected are As 
and either or (A% - Aq n ^.x)> where T) = n/n s . The 

governing equations are uncoupled and the dependent variables are 
solved one at a time at any body station in the following 
sequence: Rarikine-Hugoniot equations to obtain the shock condi- 

tions, material response for T w and m w , elemental equations 
C^, radiation transport for Q, energy equation for H, equilib- 
rium chemistry for Cp and T, s- momentum for u, global contin- 
uity for n s and v, n-momentum for p, and the equation of 
state for p. The solution is iterated until convergence is 
achieved. The solution advances to the next body station and 
uses the previous converged solution profiles as initial values 
for starting the solution at the new body station. This proce- 
dure is repeated until a solution pass is obtained. In the first 
solution pass, it is necessary to make some approximations (see 
Ref. [3]). These approximations are then removed by iterating 
the global solution. Two solution passes are generally suffi- 
cient. This solution procedure is programed for the CDC 6600 
computer. 
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RESULTS AND DISCUSSION 


To demonstrate the essential features of the analysis, 
results are presented for the following earth entry conditions: 
velocity = 15. 24 km/s and altitude = 62.2 km. The body is a 
hyperboloid with a total included angle of 120° and a nose radius 
of 0.305 ni. For these conditions, the free-stre&m Reynolds num- 
ber is 67,850 and the stagnation shock pressure and temperature 
are 0. 5 atm and 14,700 K, respectively. Mass injection rates, m, 
from 0. to 0. 60 are considered. Results for both air and 
phenolic- nylon injection are presented. For air, seven chemical 
species are used: 0, Op, 0*, N, Np, N*, and e _ . For phenolic- 

nylon, 20 chemical species are used; the seven used for air plus 
0, Cp, C 5 , CO, CN, CpH, C^H, C^H, CpHp, C + , H, Hp, and HCN. The 
present analysis can be easily modified for planetary atmospheres 
other than earth by including thermodynamic, transport, and radia- 
tive property data for the necessary species that are not pres- 
ently accounted for. 



Figure 2. Stagnation temperature profiles 
for various injection rates. 

Figure 2 shows the effect of mass injection on stagnation 
temperature profiles. Roth coupled and arbitrary injection are 
considered. The injection rate of 0.12 is the coupled mass 
injection rate for a phenolic- nylon ablator. At this point in 
the trajectory, the mass injection is sufficient to cause a 98$ 
reduction in q* and a reduction in q£. For phenolic-nylon 
injection, the specific heat of the gas is much greater near the 
wall and the radiation flux toward the wall is attenuated much 
more than for the corresponding air injection rate. Consequently, 
the phenolic- nylon injection has a much larger effect on tempera- 
ture profile and surface heating than air injection, as shown in 
Figure 2. For air injection of 0.12, the heating reduction is 92$ 
tor q* and 1 (>% for q*. For phenolic- nylon .injection rates in 
excess ot the coupled injection value, only small additional 
reductions in heating occur. 
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For large injection rates, the problem of numerical insta- 
bilities is characteristic of most, if not all, analyses that 
have been developed. For example, the analysis described in 
Reference [ 9 ] was unable to obtain a converged solution to the 
energy equation for an injection rate of 0.2. However, with the 
present analysis, a converged solution is obtained for injection 
rates substantially in excess of 0.2, as shown in Figure 2. 

Figure 5 presents the surface heating rate and mass injec- 
tion distributions where the mass injection from a phenolic- nylon 
ablator is coupled with the radiating flow-field solution. For 
the body shape considered and because the shock thickens with 
distance downstream, the radiative heating and mass injection do 
not continue to decrease with increasing distance downstream. 
Figure ^4- shows the tangential velocity profiles for the stagna- 
tion and downstream locations. These results show very clearly 



Figure 3* Shock- layer results with Figure k. Tangential velocity 
coupled ablation injection. profiles for coupled abla- 

tion injection. 

the nonsimilar behavior of the flow field. Also, note that the 
viscous effects along the stagnation streamline are evident 
throughout the shock layer. Downstream, however, the viscous 
effects are confined to a small segment of the shock- layer 
thickness. 

The spectral details of the radiation flux at the wall are 
shown for the continuum and line radiation in Figure 5 . A com- 
parison of the surface spectral flux with and without injection 
are made. The results show that for both continuum and line 
radiation, the phenolic-nylon injection causes a large attenua- 
tion of the ultraviolet radiation. 
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PHOTON ENERGY, eV PHOTON ENERGY, eV 

(a) Continuum. (b) Line. 

Figure 5. Comparison of wall radiation with and 
without mass injection. 
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